function dx=Plant(t,x,flag,para)
dx=zeros(2,1);

ut=para(1);
dt=para(2);
d=para(3);


dx(1)=x(2);
dx(2)=-0.1*x(2)-x(1)^3+dt+ut+d;